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Abstract 

We prove Euler-Lagrange and natural boundary necessary optimality conditions for fractional problems of the 
calculus of variations which are given by a composition of functionals. Our approach uses the recent notions of Riemann- 
Liouville fractional derivatives and integrals in the sense of Jumarie. As an application, we get optimality conditions 
for the product and the quotient of fractional variational functionals. 
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1 Introduction 

The present work is dedicated to the study of general (non-classical) fractional problems of calculus of 
variations. As a particular case, when a—t 1, one gets the generalized calculus of variations [1] with functionals 
of the form 



H 



f(t,x{t),x'(t))dt 



where / has n components and H has n independent variables. Problems of calculus of variations as these 
appear in practical applications (see [H-0] and the references given therein) but cannot be solved using the 
classical theory. Therefore, an extension of this theory is needed. 

The fractional calculus of variations started in 1996 with the work of Riewe [4j. Riewe formulated the 
problem of the calculus of variations with fractional derivatives and obtained the respective Euler-Lagrange 
equations, combining both conservative and nonconservative cases. Nowadays the fractional calculus of varia- 
tions is a subject under strong research. Different definitions for fractional derivatives and integrals are used, 
depending on the purpose under study. Investigations cover problems depending on Riemann-Liouville frac- 
tional derivatives (see, e. g., ta -lsj), the Caputo fractional derivative (see, e.g., |9 j - lllO . the symmetric fractional 



derivative (see, e.g., 12|, |l3|), the Jumarie fractional derivative (see, e.g. . Il4|-|2li). and others [22j-|26j. For 
applications of the fractional calculus of variations we refer the reader to 12, 14, lj| 25, 27-31|. Here we use 
the fractional calculus proposed by Jumarie. This modified Riemann-Liouville calculus has shown recently to 
be very useful in the fractional calculus of variations for multiple integrals 14], and provides an efficient tool 
to solve fractional differential equations 2l|. 

The paper is organized as follows. In Section [2] we present some preliminaries on the fractional calcu- 
lus proposed by Jumarie. Our results are then given in Section [3] We begin Section [3] by formulating the 
general (non-classical) fractional problem of calculus of variations ([3J. The problem is defined via the frac- 
tional derivative and the fractional integral in the sense of Jumarie. We obtain Euler-Lagrange equations and 
natural boundary conditions for the general problem (Theorem 13. 2p , which are then applied to the product 
(Corollary 13 A\ and the quotient (Corollary I3.6|) . In Section [4] we provide an example illustrating our results. 
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2 Preliminaries 

For an introduction to the classical fractional calculus we refer the reader to I2-34|. In this section we 
briefly review the main notions and results from the recent fractional calculus proposed by Jumarie 

Definition 2.1. Let / : [a,b] — > R be a continuous (but not necessarily differentiable) function. The Jumarie 
fractional derivative of / is defined by 

Z (q) W:=TtAt r(*-T)" a_1 (/W-/(o))dr, a<0, 



r(-«) Jo 

where F(z) = t z ~ 1 e~ t dt. For a positive a, 
in the case < a < 1, and 

/ (a) (i):=(/ Ca - n) (i)) Cn) , n<a<n+l, n>l. 
The Jumarie fractional derivative has the following properties: 

(i) The crth derivative of a constant is zero. 

(ii) If < a < 1, then the Laplace transform of f^ is given by 

£{f (a) (t)} = s a £{f(t)} -s"-7(o). 

(iii) (g(t)f(t))^= 9 ^(t)f(t)+g(t)f^(t), 0<a<l. 

Example 2.2. Let f(t)=t">, 7>0, and 0<a<l. Then f (a) {x) = T( 7 + l)r- 1 (7 + 1 - a)^" Q . 
Example 2.3. Let c and x be given constants. The solution of the fractional differential equation 

x {a) {t) = c, x(0) = x , 

is given by 

x(t) = ^-t a +x , 

where a\ := T(l + a). 

The integral with respect to (di) a is defined as the solution of the fractional differential equation 

dy = f(x)(dx) a , x>0, y(0) = y o , 0<a<l. (1) 
Such solution is provided by the following result: 

Lemma 2.4. Let f denote a continuous function. The solution of the equation fl} is 

f(r)(dr) a :=a [ (t-r)"- 1 f(r)dr, 0<a<l. 
Jo 

Example 2.5. Let /(*) = 1, and < a < 1. Then, J '(rfr) a = t a . 

Example 2.6. The solution of the fractional differential equation x^it) = f(t), x(0) =x , is 



x(t) = Xo + T- 1 (a) {t-TT- x f{r)dr. 
Jo 

In the discussion to follow, we will need the following formula of integration by parts: 

["uW(t)v(t)(dtr=a![u(t)v(t)] b a - [ b u(t)v<«\t)(dty, (2) 

J a J a 

where al := r(l + q). 
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3 Main Results 

The general (non-classical) problem of the fractional calculus of variations under our consideration consists 
of extremizing (i.e., minimizing or maximizing) 



C[x]=H^£ h{t,x{t),x^\t)){dt) a \...^ f n ((t,x(t),x^(t))(dt) a ^ 
(x(a)=x a ) (x(b) = x b ) 

over all x € D with 

D:={igC°: x^ ai \i — 1, . . . ,n, exists and is continuous on the interval[a, b}}. 

Using parentheses around the end-point conditions means that these conditions may or may not be present. 
We assume that: 

(i) the function H : W 1 — > K has continuous partial derivatives with respect to its arguments and we denote 
them by H[, i = l,...,n; 

(ii) functions (t,y,v) —> fi(t,y,v) from [a,b] x I 2 to I, i = l,...,n, have partial continuous derivatives with 
respect to y,v for all t € [a,b] and we denote them by f iy , f iv ; 

(hi) fi, i = 1,. . . , n, and their partial derivatives are continuous in t for all i£D. 

A function x G D is said to be an admissible function provided that it satisfies the end-points conditions (if 
any is given). The following norm in D is considered: 



x\\ = max \x(t)\ + V max \x (a ^(t)\. 

t£\a,b] — 't£\a,b] 



i=l 



Definition 3.1. An admissible function x is said to be a weak local minimizer (resp. weak local maximizer) 
for Q if there exists 5 > such that C[x] <C[x] (resp. C[x] >£[x]) for all admissible x with — x\\ <S. 

For simplicity of notation we introduce the operator (x);, i = l,...,n, defined by 

(x) z (t) = (t,x(t),x^\t)). 

Then, 

c[x\ = h( [ h(x),{t)(dty\..., f f n (x) n (t)(dty 



The next theorem gives necessary optimality conditions for problem ([3]). 
Theorem 3.2. If x is a weak local solution to problem ([3]), then the Euler-Lagrange equation 

n 

d H[ ( [x] , . . . , T n \x ] ) (b - 1 ) «< - 1 (/ i9 (x ) < (t) - ft ) (x ) i (t )) = 

!=1 

holds for all t € [a, b) , where T^x]^ fi{x)i(t)(dt) ai , i = l,...,n. Moreover, if x(a) is not specified, then 

n 

^a i \H^ 1 [x],...,^ n [x])f iv (x) i {a) = ] (4) 

if x(b) is not specified, then 

n 

J2 ! H[ ( T x [x] , . . . , T n [x] ) f iv (x ) 4 (6) - . (5) 
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Proof. Suppose that C[x] has a weak local extremum at x. For an admissible variation h G D we define a 
function <\> : M — >• R by </>(e) = £[(ai + e/i)]. We do not require h(a) = or h(b) = in case x(a) or x(b), 
respectively, is free (it is possible that both are free). A necessary condition for x to be an extremizer for C[x] 
is given by ^>'(£r)| £= o =0. Using the chain rule to obtain the derivative of a composed function, we get 



<l> l {e)\ e=0 = Y d H' i {J : i[x],...,T n [x]) / [/ i9 <x) i (t)M*) + /»{5> i (i)ft (ai) W] 
<=i " a 

Integration by parts (see equation @) of the second term of the integrands, gives 

f iv (xUt)h {ai) (t)(dtr* = [aJf^ixU^hmlZl- f flT\i)^t)h{t){dtr. 

J a 

The necessary condition cj>'(e)\ £=0 — can be written as 

o=53^(^i[x],...,^ n [x]) / (f iy (i)dt)-fi ai \xut))h(t)(dtr< 

i=l Ja 

n 

+ Y, H[ ( Ji [£},..., T n [x] ) [a t !/«„ (x) i (t) h (t)] 

Taking into account Lemma |2.4| we have 
•' <=i 

+ 5] [x],.. .,T n [x\) [a l lf lv (x) t (t)h(t)] t t Z b a • (6) 

In particular, equation ([5]) holds for all variations which are zero at both ends. For all such h J s, the second 
term in © is zero and by the Dubois- Reymond Lemma (see, e.g., [35j]), we have that 

n 

J2 <*iHt (Ji [£],..., T n [x] ) (b - tr - 1 (f iy (Z) i (*) - ft > (Z) i (t)) = (7) 

j=i 

holds for all £ € [a, 6). Equation © must be satisfied for all admissible values of h(a) and h(b). Consequently, 
equations ^ and ([7]) imply that 

n n 

= ^J^ 1 [s],...,j;[5]te^ (8) 

1=1 i=i 

If x is not preassigned at either end-point, then h{a) and h(b) are both completely arbitrary and we conclude 
that their coefficients in © must each vanish. It follows that condition ^ holds when a; (a) is not given, and 
condition ([5]) holds when x(b) is not given. □ 



Note that in the limit, when a t — > 1, i = l,...,n, Theorem 13.21 implies the following result: 
Corollary 3.3 (Th. 3.1 and Eq. (4.1) in [1]). If x is a solution to problem 



£[x}=H^J / 1 (t,a;(t),a;'(t))^,...,y f n (t,x(t),x'(t))dt\ — > extr, 
(x(a) = x a ) (x(b)=x b ) 

then the Euler-Lagrange equation 

H I {T x [x] , . . . , T n [x] ) U v (t, x (t) , x' (t)) - -^-f iv (t, x (t) , x'(t))\=0 

1 = 1 v ' 



A.B. Malinowska, M.R. Sidi Ararai and D.F.M. Torres 5 

holds for all t £ [a, b], where Fi[x} = J b fi(t,x(t),x'(t))dt, i = l,...,n. Moreover, if x(a) is not specified, then 

n 

^2 H 'i(^ [x],...,jF n [x])f iv (a, x(a) , x' (a)) = 0; 

i=l 

if x{b) is not specified, then 

n 

]T H[(F 1 [x],...,F n [x])f lv (b, x(b) ,x'(b)) = 0. 
1=1 

Corollary 3.4. If x is a solution to problem 

C[X]= (/ 6/l(i)lW(cft)ai ) (/ M^2(t)(dt) a A ^ extr, 
(x(a) = x a ) (x(b)=x b ) 

then the Euler- Lagrange equation 

holds for all t£ [a,b). Moreover, if x{a) is not specified, then 

oil ! Ti [x] fiv (x ) i (a) + a 2 \ T x [x] f 2v (x) 2 (a) = 0; 

if x(b) is not specified, then 

ai ! T 2 [x] fi v {x) i (6) + a 2 ! Ji [x] /20 (x) 2 (6) = 0. 

Remark 3.5. In the case — > 1, i = 1,2, Corollary 13.41 gives a result of [lj: the Euler-Lagrange equation 
associated with the product functional 



£ M=(j h{t,x{t),x'{t))dtj[J f 2 {t,x(t),x'(t))dt 

is 

F 2 [x] (^f ly (t,x(t),x'(t))- ^f lv (t,x(t),x'{t))^ +F 1 [x] ^f 2y (t,x(t),x'(t))-j t f 2v (t,x(t),x'{t)^j =0 

and the natural condition at t = a, when x(a) is free, becomes 

T 2 [x] f lv (a, x (a) , x' (a)) + T x [x] f 2v (a,x(a) ,x'(a)) = 0. 
Corollary 3.6. If x is a solution to problem 

£[*] = \ > extr > 

f b J 2 (x) 2 (t)(dt)<** 
(x(a) = x a ) (x(b)=x b ) 

then the Euler-Lagrange equation 

a^b-tY-- 1 (hyixUV-ftHzhit)) -a 2 Q(b-tr^ (f 2y (x) 2 (t) - '<2>a(f)) = 

holds for all tE [a, b), where Q — . Moreover, if x{a) is not specified, thena 1 \f lv {x) 1 {a)—a 2 \Qf 2v (x) 2 {a)={); 
if x(b) is not specified, then a 1 \f lv (x) 1 (b) — a 2 \Qf 2v (x) 2 (b)=Q. 

Remark 3.7. In the case a t — > 1, i = 1,2, Corollary 13.61 gives the following result of [l|: the Euler-Lagrange 
equation associated with the quotient functional 

= tfh(t,x(t),x'(t))dt 
J b a f 2 (t,x(t),x'(t))dt 

is 

f lv {t,x{t),x\t))-Qf 2v {t,x{t),x\t))--[{f lv {^ 
and the natural condition at t = a, when x{a) is free, becomes 

f lv (a,x(a),x' (a)) - Qf 2v (a,x(a),x' (a)) = 0. 
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4 An Example 

Consider the problem 

minimize C[x]=(^J {x^\t)f{dt)^]( j t?x^\t){dt)? 



(9) 



x(0)=0, x{l) = l. 

If x is a local minimizer to ([9]), then the fractional Euler-Lagrange equation must hold, i.e. 

ig 2 (l-i)-52(5(i'(0) ( * ) + ig 1 (l-t)^(^) ( i ) =0, 

where i 

g 1= / (z*(t)) 2 (dt)3, Q 2 = / tix(i'(t)(dt)l. 



Hence, 

Q 2 2(^i)(t))^ + Q 1 ^ = 0. (10) 
If Q 2 = 0, then also Qi = 0. This contradicts the fact that a global minimizer to the problem 

minimize F x [x]= I (x ( ^ (t)) 2 (dt)i 
Jo 

x(0)=0, x(l) = l 

is x(t) = th and F^x] = (^) 2 . This can be easily shown by the results obtained in [15| . Hence, Q 2 ^ and 
(flU)) implies that candidate solutions to problem © are those satisfying the fractional differential equation 

(*(§)(*))(*) (n) 

4^2 

subject to the boundary conditions x(0) = and x(l) = 1. Solving equation (jTTJ) we obtain 

X ^^V^L V 8Q 2 2Q 2 

Substituting (|12p into functionals and J-" 2 gives 

f 1_ -32Q! 2 -4n .Q^tt 

= Q 2 

We obtain the candidate minimizer to problem Q solving the system of equations (fl 

l - 85 — ^ r2 J (< - r) 2dT ' 

where 



i rtQ^+w* % r i] (i - T) -f d , (12) 



2 2 

Q 2 2 11 = < 5 ] 

J_ -32Qi+3Qi7r 2 + 12 7r 3 / 2 C 
. 96 Q 2 



(13) 



4 7r(0F(i7ri + i V^ 3 - 8tt) - 4) 



Q » = 3 -32^ ^ °' = 12^+12^ 
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